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Abstract: We study the Gaussian hermitian random matrix ensemble with an 
external matrix which has an arbitrary number of eigenvalues with arbitrary 
multiplicity. We compute the limiting eigenvalues correlations when the size of the 
matrix goes to infinity in non-critical regimes. We show that they exhibit universal 
behavior and can be expressed with the Sine and Airy kernels. We also briefly review 
the implication of these results in terms of non-intersecting Brownian motions. 
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1 Introduction and statement results 

In these notes, we are interested in the hermitian matrix model in an external field 
defined by the partition function: 

Z:= f dMe~ NT <^- AM ) (1-1) 

J H N 

where one integrates over the hermitian matrices M of size TV x N, A is a deterministic 
diagonal matrix of the form 

ni ri2 n k 

A := diag(ai, . . . , a u a 2 , . . . , a 2 , . . . , a k , . . . , a k ) (1-2) 

with di 7^ ctj for i ^ j and the measure dM is the product of the Lebesgue measure of 
the real entries of the matrix M: 

dM := ]JdRe (Mij) dim (A%) \\dMu. (1-3) 

i<j i 

The main motivation of this work comes from the study of the exclusion process de- 
fined as follows. Consider iV independent non- intersecting Brownian motions {^(t)}^ 
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starting from x = at time t = and ending at point x — Oi(l) at t — 1 for of them 
with z = 1, . . . , k. One obviously imposes 

Ei=i n» = N. Indeed, it is known [TQl d 

that, up to a rescaling, 

mi = ^ \ ? (1-4) 

the position of the particles £,(£) at any time t G (0, 1), are distributed as the eigenval- 
ues vrii of the random matrix M submitted to an external matrix A(t) whose eigenvalues 
are given by 

Oi(*) = Oi(i)y^. (i-s) 

We study the behavior of this process when the number N of Brownian motions gets 
very large, i.e. the size of the random matrix M gets large and the filling fractions 
= 6i are kept fixecjE In this limit, at a given time t the Brownian motions accumulate 
to fill a set of l(t) intervals with 1 < l(t) < k and the behavior near the boundaries: 
/(£) — > 1 as t —>■ and l(t) — > k as t — > 1. Indeed, at the beginning, all the particles 
start at the origin and stay next to x = at short times. At the end, the particle are 
split into k groups, each composed of n, particles, lying next to the k and points aj(l). 
As time grows, the unique interval corresponding to t = splits into two intervals for 
a critical value of time corresponding to a phase transition. Keeping time growing, 
these intervals split also for critical values of time and this process keeps on iteratively 
until one gets k intervals reaching the end points. In these notes, one describes the 
distribution of Brownian particles in these intervals for any non critical value of time by 
studying the associated matrix model for different regimes. In section [3j we study the 
case of large time when l(t) = k as a warm up before studying the generic non-critical 
case 1 < lit) < k in section HI 

In both section, we follow the same procedure. We first study the Pastur equation 
[9], or spectral curve, 

^_ a , (i)) ^_,_|:_^_) =0 , (1 _ 6) 

of the matrix model obtained for time t, as well as the associated Riemann surface. 
The structure of the Riemann surface changes with time: especially, it has as many 
cuts as the number / of intervals [z 2 i-i, z 2 i\ i=l filled by the eigenvalues. We then study 
an associated Riemann-Hilbert problem of size (k + 1) x (k + 1) whose solution Y(x) 



2 In order to impose this condition, one has to restrict our study to ./V common multiple of the 
denominators of the fractions . It should be possible to release this constraint by showing that the 
corrections one should get by adding an isolated eigenvalue to A are not seen at the leading order 
of the physical quantities when N — > oo. This issue is not essential for the present work and is not 
studied here. 
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allows to define a fundamental kernel Kn(x, y) giving access to any correlation function 
of the eigenvalues through: 

R m (x 1 ,X2,...,x m ) = <let(K N (x j ,x k )) 1 ^ jk ^ m . (1-7) 

By successive transformations of this RH problem, we turn it to a simple one expo- 
nentially close when iV — > oo to a model RH problem whose solution can be written 
explicitly. This asymptotic solution of the RH problem leads to the main theorems of 
this paper. The first one describe the limiting density of Brownian motion at a given 
time when we have a large number of them. 

Theorem 1.1 The limiting density of eigenvalues 

K N (x,x) 



p(x) := lim 



(1-8) 



exists, is supported by the real cuts [z2i-i, z 2i\\ =1 and is expressed in terms of the unique 
solution yo(x) of the Pastur's equation 



t[(y-<k(t)) 



x-y 



E 



n j 



i=i 



-N(y- aj (t)) 



= 



such that 



by 



y (x) ~ - + 0(\) 



as 



x — > oo 



x 



X" 



1 



(1-9) 

(1-10) 
(1-11) 



p(x) = - \lmy (x) \ . 

TV 

It is analytic on the cuts and vanishes as a square root at the edges Zi. 

The second one describes the behavior of the kernel inside the support of the eigenvalues 
(or Brownian motions). 

i 

Theorem 1.2 For any x E [_J(z 2 i-i, z 2 i) and any pair (u,v) E M 2 ; the kernel converges 

i=i 

to the universal sine kernel: 



lim — - 

N^oo Np(x) 



K N [x + 



u 



-, x + 



Np(xY Np(x) 



sin tt(u — v) 
tv(u — v) 



where we use the rescaled kernel: 

K N (x,y)=e N ^- h ^K N (x,y) 

with the function 



x 



h(x) := -— + Re 
where the primitive will be made precise later. 



yo(x) 



(1-12) 

(1-13) 
(1-14) 
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The last one describes the asymptotic behavior of the kernel at the edge of the support 
of eigenvalues Zj. 

Theorem 1.3 For any i = 1, ... ,2k and every (u, v) EM 2 , the rescaled kernel is given 
by the universal Airy kernel on the edges of the supports of eigenvalues: 

km T K N Zi + T ,zi + j = — — — 

(1-15) 

where pi is a constant characterizing the behavior of p near the branch point z^: 

p( x ) = P±\ x - Zi \k(l + {l)) (1-16) 

7T 

when x Zi. 

Organization of the paper: 

• section [2] is a reminder on the Gaussian matrix model with an external field and 
its link to non-intersecting Brownian motions. 

• section [3] is a warm up where we study the asymptotics of the kernel and mean 
density of eigenvalues in a particular non degenerate case: there are as many 
supports of eigenvalues as the number of different eigenvalues of the external 
matrix A: I = k. 

• in section HJ we study the general case 1 < I < k and prove the theorems II. H 11.21 
and 11.31 for any non-critical spectral curve Eq. fll-6p . 

• section [5] is a conclusion and present some possible generalizations of the present 
work. 

2 Gaussian matrix integral in an external field, re- 
minder 

2.1 Definitions 

In this section we consider the hermitian matrix model defined through the partition 
function 

Z:= ! dMe~ NT <^- AM ) (2-1) 

Jh n 

where one integrates over the hermitian matrices M of size iV x N, A is a deterministic 
diagonal matrix of the form 

A := diag(ai , . . . , ai, a 2 , ■ ■ ■ , a 2 , . . . , a k , . .. , a k ) (2-2) 
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and the measure dM is the product of the Lebesgue measure of the real entries of the 
matrix M: 

dM := Y[dRe (A%) dim (A%) JJdM a . (2-3) 

i<j i 

By diagonalizing the matrix M, we can reduce this problem to the study of its 
eigenvalues x%, . . . ,x^, which are constrained by the probability measure: 

dp(x u ...,x N ) = Y[ dxi^re- N ^ i det (e"^) (2-4) 

where A(x) denotes the Vandermonde determinant and one considers the eigenvalues 
di with their multiplicities. 

They can be seen as a Coulomb gas composed k groups of particles living on 
the real axis and submitted to a gaussian potential Vi(x) = ^ — aix and a logarithmic 
repulsive interaction In \x% — Xj\. The position of the eigenvalues is then an equilibrium 
configuration for such particles. When the size of the matrix tends to infinity, the 
number of eigenvalues grows in the same way as the latter and they fill a finite set 
of intervals [z2i-i, Z2if i=1 with 1 < I < k. The number I of such intervals as well as 
their positions depend on the a^s and the ratios 6j := ^. They can be encoded in an 
algebraic curve as we remind now. 

Let us define the density of states as the one point correlation function: 
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N 



p(x) := -j- / d{i(ui,u 2 , . . . ,u N )y^5(x - Ui), (2-5) 
JV Jr n ~[ 

as well as its Stieltjes transform, called the resolvent 

W{x) := [ ^-du. (2-6) 
Jr x ~ u 

Remark that the density of states p is expected, by the analogy with a Coulomb 
gas, to have a compact support composed of the / intervals corresponding to the / 
saddles of the effective potential for the Coulomb particles. Thus, the resolvent can be 
written: 

W(x) = J2±-[ J^du (2-7) 

and is continuous except cuts on the support Ui[z 2 i-i, z^] of the density p. 
One can show [SI El [7] that this resolvent satisfies an algebraic equation: 

S(x,W(x)) = (2-8) 
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where S(x, y) is a polynomial in its variables. Once the so-called spectral curve 
S(x, y) = is known, one thus obtain W(x) by solving a polynomial equation. Never- 
theless, generically there exists as many solutions to this equation as the degree of this 
polynomial in y. But, for some particular values of x, this equation has double roots 
in y. These particular points Zi are usually called branch points since they correspond 
to the points where two different solution of the equation meet. 

It is easily seen that these branch points correspond to the Zi fixing the endpoints 
of the different component of the support of the density p. Finding the position of 
these branch points allows thus to know the support of the eigenvalues. 

When the size of the matrices tends to infinity, the eigenvalues fill this support 
continuously and one is interested in the behavior of the level density inside these 
support and on their boundaries. This is why we now study this spectral curve in more 
details. 

2.2 Algebraic curve 

Following |7J, one can compute the spectral curve associated to the considered model. 
For this model and arbitrary values of the a^'s: the algebraic equation takes the form: 

k 

£(x, y) := (x - y) \[{y - a t ) + P(y) = (2-9) 
i=i 

where P{y) is a polynomial of degree k — 1 whose coefficients must be fixed by the 
behavior of the function y at infinity. 

Indeed, this algebraic equation can be seen as the realization of a compact Riemann 
surface £ embedded in CP 1 x CP 1 with two meromorphic functions x and y such that 

VpGS, S(x(p),y(p)) = 0. (2-10) 

One can easily see that this is a genus curve composed of k + 1-sheets in y. It 
means that for a given value of x(p) E CP 1 , there exist k + 1 points on S which are 
pre-images of x. We note them p ^ 1 , . . . ,p k : 

Wi = 0,...,k, x(p i ) =x(p). (2-11) 

To agree with the notations of [3j HJ E], one can define the value of the function y 
on the different sheets as k + 1 different functions on the complex plane of the image 
x{p): 

Definition 2.1 Let the functions : C — ► C be defined by 

Ux(p))--=y(p 1 )- (2-12) 
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The different sheets can be discriminated by the behavior of y(p) when x(p) — * oo. 
From the definition, one can label the different sheets in such a way that: 

y(p°) ~*(p)->oo x{p) + ^ + O (^J 
Vz + , yip') ^(pj-.oo a, + f + O 

i.e. the ^-functions have the following asymptotics: 



x(p) 2 



(2-13) 



^0(X)~™X+ J +0 (4y) 



(2-14) 



Inserting these asymptotics inside the algebraic equation Eq. (12—91) . one gets the 
coefficients of the polynomial P(y) and the algebraic curve reads: 




(2-15) 



2.3 Kernel and correlation functions 

In [HI E], Daems and Kuijlaars extended classical results of random matrix theory to the 
case of 1 matrix model in an external field: they explain how the correlation functions 
can be computed thanks to a unique kernel which can be linked to the solution of an 
associate Riemann-Hilbert problem. 

They showed that, using some multiple orthogonal polynomials, one can express the 
correlation functions of the model (see [SI IB] [HI H2] ) • In particular, the joint probability 
density on the eigenvalues can be written 

^ydet(KAr(x i ,x fc ))i< ijfc <iv (2-16) 
and the m-point correlation functions 

Rmixi,..., x m ) = det(K N (xj,x k )) 

l<j,k<m 

(2-17) 

where the kernel K^(x,y) is expressed through some multi-orthogonal polynomials if) 
and by 

JV-l 

K N (x,y) = e-^ x2+ ^J2Mx)Mv)- (2-18) 

fc=0 

Moreover, in [SI IB], it was shown that these kernels associated to multiple orthogonal 
polynomials are solution of a general Riemann Hilbert problem. Let us specify here 
the problem whose solution is the kernel giving correlation functions of the hermitian 
matrix model introduced above with an external matrix A whose eigenvalues are Oj, 
i = 1, . . . , k. 
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Theorem 2.1 (Daems-Kuijlaars) 

The correlation functions are given by 

R m (xi, . . . , x m ) = det(K N (xj,x k ))i<j,k<r 
where the kernel Kjsr(x,y) reads: 



K N (x, y) 



n x {y)Y7\y)Y + {x)^{x) 



2in(x — y) 

where we define the vectors of size k + 1 : 

n^x) = [0 ^^(j) . . . e N ^u k {x)] 
with Ui := e - N (\~ a i x ) an d 

Q 2 {x) = [1 ... 0], 
and Y(x) is the unique solution of the following Riemann Hilbert problem: 

• Y : C \ R -> C( fc+1 ) x ( fe+1 ) is analytic; 

• for x 6 M, one has the jumps 



Y+(x) = Y-(x) 



( 1 Ui(x) UJ 2 (x) 

1 
1 



\ 







u k (x) \ 



1 J 



where Y + (x) and Y-(x) denote respectively the limit ofY(z) when z 
the upper and lower half planes; 



• when x — > oo ; one has the asymptotic behavior: 



Y(x)=[I + 0[- 



( z N 

z~ n ^ 
z'^ 



o \ 






\ 

where I denotes the identity matrix of size (k + 1) x (k + 1). 



/ 



2.4 Non-intersecting Brownian motions 

One of the main motivations for this work is the study of a physical system composed 
of many N Brownian motions as N gets large. Let us consider N Brownian motions 
Xi(t), i = 1, . . . , N, on the real axis starting from the origin at time t — 0: 

Vi = l,...,7V, Xi (0) = 0. (2-25) 

Let us impose that they cannot intersect and that their final positions at time t is given 
by a set of k ordered distinct real points af 

(ni ni rife \ 

a77"^- , ai, a 2 , . , ai, . . . , a k , . ?. , a k ) (2-26) 

for any integers {n,i} k i=l such that 

k 

^m = N. (2-27) 

i=i 

The physical quantities describing this statistical system at a given time < t < 1 are 
the correlation functions: 

R m (X 1 ,X 2 ,...,X M ;t) := — - / dP( Xl (t),x 2 (t), x N (t)) TT - X^) 

(2-28) 

where P(xi(t), x 2 (t), . . . , x^it)) is the probability measure on the Brownian motions 
at time t. 

This problem does not look easy to solve at first sight. Fortunately, there exists 
an underlying integrable structure and it was proved that these correlation functions 
coincide with the correlation functions of a hermitian gaussian matrix model submitted 
to an external matrix 

ni ri2 rife 

/V A. /V 



A := A(t) := diag(ai(t), . . . , a^t), a 2 (t), a 2 (t), a k (t), a k (t)) (2-29) 
where the eigenvalues dj are just rescaled versions of the ending points: 



u,(t) :=, M /^1_. (2-30) 

In order to compute these correlation functions at time t, one is then left with the study 
of the corresponding random matrix M(t) whose statistic depends on t only through 
the external matrix eigenvalues a«(t). Remark also that the time evolution of these 
external eigenvalues is very simple. 
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Let us now interpret the ingredients of the matrix model's study in this statistical 
physics setup to get some intuition. First of all, one sees that the eigenvalues of the 
random matrix X{(t) are identified with the position of the Brownian motions at time t. 
The cuts of the spectral curve [z2i-i{t), Z2i{t)] l i=1 should then coincide with the support 
of the Brownian motions at time t. One can then study the evolution of these support 
by studying the time evolution of their edges, i.e. the time evolution of the branch 
points. 



Consider the equation o 
giving x in terms of y read 



the algebraic curve Eq. (I2-I5p . The inversion formula 




(2-32) 



and has, generically, 2k branch points given by the roots of the equation: 

k 



1 = Y — 3 , ^ . (2-33) 

Let us call {^i(i)}|=i these branch points ordered following their real part. Since the Oj's 
are real, the roots of Eq. (12-331) are either real or come in pairs of complex conjugated 
points. 

Let us now have a look at the two extremal cases corresponding to the boundaries. 
For t — > 1, the eigenvalues diverge a«(t) — > oo. Thus, the equation for the branch points 
Eq. (12-331) as solutions: 

y ~ Oi(l) ~ oo, (2-34) 

i.e. one has 2k real branch points of the type Zi%-\ < ctj < Z2% with Z2i~i and Z2% 
very close to a^. This confirms the physical intuition of the behavior of the Brownian 
motions: they merge in k groups close to the end points dj when t — > 1. On the 
contrary, as t — > 0, all the a,i(t) merge to and all the branch points merge also to 0: 
only two of them are real, the other 2k — 2 come in pairs of complex conjugated points. 

Let us now come to the generic configuration. Starting from t — 1, as the time 
decreases the real branch points Z2i-i and Z2i move away from a^. When Z2% and Z2«+i 
meet, one reaches a critical point and the corresponding configuration is not studied 



^Remark that this gives a natural rational parametrization of this curve: 

f V = z 



k ru (2-31) 



X - Z + J2i=l N{z- ai (t)) 
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in this paper but will be addressed elsewhere. After crossing this phase transition, 
these two branch points move away from the real axis becoming complex conjugated. 
Then, at any intermediate time, a generic configuration consists in 21 real branch points 
(^2i-2) z 2i), i — 1) • • • > I, where the number of real cuts / decreases with time from k to 
lc. The other branch points are complex conjugated numbers {wf'\w < f 1 ^ whose real 
part belong to the segments [z 2 i-2, Z2i\- 

In the following sections, we study the Riemann surface associated to these different 
configuration: in section El we consider the case where we have only real branch points 
whereas in section HI we study the generic case with an arbitrary number of real 
eigenvalues. 



3 Large time case 

In this section, we assume that all the branch points are real. This configuration is 
encountered for large times in the brownian motion setup: it corresponds to having all 
the particles in the neighborhood of their end points. Then one can simply mimic the 
computations of Bleher and Kuijlaars [5] and study the algebraic curve. 

3.1 Study of the spectral curve 

Before solving the Riemann-Hilbert problem, one needs to study the associated spec- 
tral curve more carefully. Remember that, in the large time case we consider in this 
paragraph, the spectral curve has 2k real branch points {zi}^. 

We can show the following lemmas concerning the sheet structure of the spectral 
curve built out of the y function. First of all, we can build analytical continuations of 
the & as follows: 

Lemma 3.1 The functions ^ for i / (resp. £o) are analytic in C\[z2i-i, z 2 i] (resp. 

k 

C \\J[ Z 2j-l,Z 2j })- 
3=1 

proof: 

It is easily proved by starting from the asymptotics Eq. ( 12- 14ft and analytically 
continuing these functions using the rationnal parametrization of this genus curve. 
□ 

These functions satisfy the following jump conditions: 

4 The number I corresponds to the number of disjoint segments filled by the Brownian motions at 
a given time t. The proof of this property is not technically involved but rather long, that is why it 
was decided not to include it in the present paper since the interested reader can easily derive it. 
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Figure 1: Sheeted structure of the spectral curve in the large time case. The sheet is 
linked to the other k sheets by non intersecting segments. 



Lemma 3.2 For any i = 1 



. . ., k and x G [z 2 i-i, z 2i ] 



£,o+(x) = Co-fa) = &-(ar) = &+(x) 



(3-1) 



and for x G UiLil^-i, z 2i [ 



Im^ +(a;) > 



(3-2) 



proof: 



It also follows directly from the definitions. □ 

It means that the spectral curve is composed of k + 1 sheets linked by k real cuts 
as follows: the sheet i is linked to the sheet by the cut [z 2 i-i, z 2 i\ (see figH]). 

We can then define the "eigenvalue density" corresponding to the density of brow- 
nian particles at a given point x of the real axis: 

Definition 3.1 Let us define, for any x G [Jj = i[z 2 j-i, z 2 j], 



p(x) := -Im£o+(aO. 

7T 

One can now study its properties giving the physically expected result: 



(3-3) 



k 



Lemma 3.3 For x G Mj^i-i, z 2 \, p(x) > and 




i=i 




(3-4) 
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On the edges of the spectrum, there exists pi > such that 

p(x) = —\x — Zi\(l + 0(z — zA) as x —>■ Zi. (3-5) 

7T 

proof: 

From the preceding lemma, one directly get that 

k 

Mx E {J}z 2i -i, z 2i {, p{x) > 0. (3-6) 

i=l 

Consider now any integer i G [l,k] and x G]z2i-i, %h [> one has 



p(x) = -Im£ 0+ (x) = —(£ 0+ (x) - £ 0+ (x)) = — (£i-(x) - £i+(x)) (3-7) 

7T ZlTl LlTX 

and thus, one can compute the integral 

P{X) = ^~\ 21 " 6+0*0) = ^~{ ^ X ) dx ( 3 " 8 ) 

where the integration contour Tj encircles the cut [z 2 i-\., z 2 j\. Since £j has no other 
poles except infinity, one can move it and use the asymptotics: 

p(x)dx = — Res £i(x)dx = — \. (3-9) 



The last part of the lemma follows from the usual behavior on the edge of the 
spectrum (see for example [3]). □ 

Remark 3.1 All the results of this lemma simply follows on the behavior of he 1-form ydx 
in the neighborhood of a cut linking two simple branch points. Thus all these local result 
could be obtained by the former study of simpler models where one has only two sheets for 
example. 

We now define the integrals of the £j functions as: 
Definition 3.2 We define the integrals Xj~ by 

Ajfc(z) := / Zk(x)dx (3-10) 



where the integration constants are constrained by: 

A (^2fc) = A fc (z 2 fc) = (3-11) 

and 

Vi = 1, . . • , k - 1 , Xi{z 2i ) = \+{z 2 i)- (3-12) 



14 



Lemma 3.4 With this definition, the jumps and asymptotics of these X-functions are 
given by: 

• When z — > oo: 



^--\nz + l + O (- 
2 V z 

and, for i ^ : 



X (z) = --\nz + l + O(-) (3-13) 



\{z) = a iZ + ^ In z + k + O ( \ ) ( 3- 1 I ) 



N — ■ -» ■ - ^2 

where U 's are some constants chosen to satisfy the constraints Eq. A3-ll\) and 
Eq. BAm . 

• They follow the jumps: 

X 0+ (x) - Ao_(ac) = -2in J2j=i+i n x e z 2i+i] 

A +(x) — A -(x) = — 2i7T X G ( — OO, Zi] 

Vi G [1, fc] , A i+ (x) - Ai_(a;) = 2i7r^ x G (-oo, z 2 i-i] 

A 0+ (x) - \ k -(x) = , X -{x) - X k+ (x) = x G [-sffflb-i, z 2 fc] 

Vi G [1, fc - 1] , A +(x) - Xi-(x) = 0, Ao-(a:) - A i+ (x) = 2iirJ2j=i+i n x E [z 2i -i, z 2i ] 

(3-15) 

proof: 

The first part of this lemma is just the rewriting of the integration of the asymptotics 
of ydx on the different sheets. 

The second part is obtained by moving the integration contours so that one is left 
by the computation of a residue at infinity on the different sheets. □ 

Let us now state a fundamental theorem for the following 

Theorem 3.1 For any i = 1, . . . , k, on R\[z 2 i-i, z 2i ], one has ReAj + < ReAo-. 
For any i = 1, . . . , k, there exists a neighborhood Ui of]z 2 i-i, z 2i [ such that 

Vj ^i, ReXj{z) < ReX (z) < ReXi(z) (3-16) 

for any z G Ui-]z 2 i-i, z 2i [. 

proof: 

Once again, one can follow |3J for the first part of the theorem: for x > z 2 i, £o(x) > 
£i(x), and Xo(x) > Xi(x) because Xi(z 2 i) = Xo(z 2 i) and A£ = ^; for x < Z2i-i, one easily 
sees that Re£ Q (x) < which gives the result. 

For the second part, one has to study the properties of the A-functions more care- 
fully. Consider the function Fi := Aj + — Ao+. It is purely imaginary on ]z 2 i_i, z 2 \ and 
its derivative F\x) = —2mp{x) has negative imaginary part. Thus ReAo(z) < ReXi(z) 
in a neighborhood of ]0 2 j_i, z 2i [. The other inequality follows from the definition and 
continuity. □ 



15 



3.2 Riemann-Hilbert analysis 

In this section, we solve the Riemann-Hilbert problem leading to the computation of 
the kernel. We look for the unique function Y defined as follows: 

Definition 3.3 Let Y be the solution of the following constraints: 

• Y : C\R^ c( fc+1 ) x ( fc+1 ) is analytic; 

• for iGl, Y(x) has the jumps 



Y:(x) = Yjx) 



( 1 uji(x) uj 2 (x) 
1 
1 



V o 











u k (x) \ 



1 J 



(3-17) 



where Ui(x) := e ~ N (—- a i x ) Y + (x) and Y^(x) denote respectively the limit of 
Y(z) when z — > x from the upper and lower half planes; 

when x — > oo ; one has the asymptotic behavior: 



Y(x) 



I + O 



/ 


z N 





. 


. 


\ 









. 


. 












z~ n2 . 


. 




\ 








. 


. z~ nk 


J 



(3-18) 



The kernel can then be written: 



K N (x,y) 



2m{x — y) 



n 1 (y)Y- 1 (y)Y + (x)n t 2 (x) 



(3-19) 



with the notations of Eq. f)2-2ip and Eq. (12-221) . To do so, we transform the problem 
step by step to turn it to a simple one (see [2] for a review on the subject). 

It is important to keep in mind the constraints we impose on the n'^s and N: 



We fix a set of k positive rational number called filling fractions: 



Pi 

e. = — 

Qi 



(3-20) 



where pi and are coprime and 



E 



i=i 



(3-21) 
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iV is a common multiple of the g^'s. 
The integers n« are fixed by the condition: 

m = eiN. 



(3-22) 



With these constraints, the spectral curve associated to this RH problem does not 
change as N grows. 



3.2.1 First transformation 

Using the notations of section 13.11 we can define 

-Nl a -Nh 



T(x) := diag(e J '"",e 



e- m >) F(x)diag ( e^^-V^^Ax-^) 



(3-23) 



N(\ k 



for x G C\R. 

Then one can check that it is solution of the Riemann-Hilbert problem: 

• T is analytic on C\R; 

• for x G R, one has the jumps 

T + (x) = T4x)j T (x) 

where 



(3-24) 



/ e Ar(A 0+ (x)-A _(a;)) g iV(A 1+ (x)-A _ (x)) e N{X 2 +(x)-X -(x)) 
q e N(\i+(x)-Xi-(x)) q 



3t(x) 













D N{X 2+ {x)-X 2 -(x)) 



when x — > oo, one has the asymptotic behavior: 



T(x) 



I-i ( )[- 

x 



e N(X k+ (x)-X -(x)) \ 




e N{X k+ (x)-\ k _(x)) 

(3-25) 



(3-26) 



Note that one can write the jump matrix Jt{x) depending on the cut on which x 



lies: 



For x G [z2i-i, Z2i], jr takes the form 



( e N(A (x)-A;(x)) + e JV(A 1+ (x)-A _ (x)) 
1 



f ,N(A i _ 1+ (x)-A _(x)) 




e «(A (x)-A i (x))_ 




e N(A;_ 1+ (x)-A _(x)) 



(3-27) 



17 



For x G R\^J[^2i-ij z 2i], 3t reduces to 



Jt 



/ i e N(Xi+(x)-X -(x)) e N{A2+(x)-\ -(x)) 

1 

1 



\ 
The kernel can then be written: 







e N(\ k+ (x)-\ -(x)) \ 






K N (x,y) 



' 4 4 



-0 ljT (y)Tr 1 (y)T + (x)0* (x) 



where 



and 



2i7r(x — y) 
= [0 e 

n 2T (x) = [e~ NXo +W o ... 0]. 



7VA 1+ (x) _ _ e NX k+ (x)i 



3.2.2 Second transformation and opening of lenses 

Remark hat one can further factorize the matrix j't(^) for x G \z%-\-, Z2i}'- 

3t(x) = ji,T(x)j s (x)j iiT (x) 
where, for x G [z 2 i-i, z 2 i], one has defined the k + 1 x k + 1 matrices: 

3s(x) = [{js)ap\ at p = Q '■= [(1 — $ai — £qo)<W + S a0 S/3i — & ai& po[ a> p = Q > 



AT(Ao-Ai)- 



<^' e 



N(Xj—Xi)- 



a,(3=0 



and 



(3-28) 



(3-29) 

(3-30) 
(3-31) 



(3-32) 

(3-33) 
(3-34) 

(3-35) 



a,/3=0 

This allows us to define the second transformation. Consider a set of k lenses with 
vertices 2:21- 1 and (see fig!2D contained in the neighborhood [/, defined in lemma |3~T1 
and define 

T(z)Ji(x) in the upper lens region 
T(z)Ji(x) in the lower lens region 



S(z) :- 



(3-36) 
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Im(x) 
▲ 




Z 2i -2 Z2i-K 



z 2i+l 



► Re(x) 



Figure 2: Lenses used in the first transformation: the lense i has vertices zn-\ and z 2i 



with 



and 



JAx] 



a,p=0 



Ji(x) : = 
for z in the lens i , and 



(3-37) 



(3-38) 



a,/3=0 



S{z) := T(z) (3-39) 

for z outside the lenses. 

We now define the jump matrix js(x) on all the discontinuity contours of S(x): 



• On [z 2 i-i,z 2i \, for % = 1, . . . , k: 

• On the upper boundary of the lens z: 

is 0*0 



(3-40) 



*afl-UE^/ ( * rA,) ) 



(3-41) 



a,P=0 



On the lower boundary of the lens i: 



(3-42) 



a,/3=0 



• OnR\\J[z 2i - 1 ,z 2i } : 

J5(x)=Jt(x). (3-43) 
Using these definition, one sees that S(x) is solution of the following RH problem: 
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• S(x) is analytic on C\(R(jr) where T denotes the boundaries of the lenses; 

• S(x) has jumps: 

Vz G R (J T , S+ (x) = S- (x) js (x) ; (3-44) 

• when x — > oo: 

5(x) =1 + f ~) (3-45) 

and the kernel can be expressed in terms of this matrix S by the following expression: 
for (x,y) G (z 2i -i,z 2 i) 

N[h(y)-h{x)} 

K N (x, y) = ^ - - (3-46) 

2 

/i(x) = ~ + Re [A 0+ (x)] , (3-47) 

i—l k—i 

Q 1Ai (x) = [-e Nilm ^+^ 0~^~0 e - m ^o+(x)] ^^o] (3_ 48 ) 



where 



and 

i— 1 fc— i 

fi2,s,i(a:) = [e-^^W] (H?~"o e mim[Xo+{x)] 0~^~0\. (3-49) 

One can see that, as N — > oo, this RH problem reduces to a model RH of the same 
kind of that studied in [3J 0J E] thanks to the properties of the A-functions in section 
13.11 Let us solve it in the same way. 

3.2.3 Model RH problem. 

This means that we want to solve the following RH problem: 



M(x):C\|J[^-i,^]^C fe+1 + 

i=l 



M(x) is analytic in C\^J[z 2 i-i, Z2i]] 

i=l 

M has jumps on the cuts [221-1 j z 2i] given by: 

k 

Vx G (Jfe-i, 2a] , M + (x) = M4x)j s (x); (3-50) 



i=l 

as x —> 00: 



M(x) =1 + 01-1. (3-51) 
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Let us lift this problem to the associated spectral curve by cutting the complex 
plane along the different y-sheets: one defines: 



for i = 1 , . . . , k and 



= ^i(C\[z 2 i-i, z 2i \) 



^0 = £o(C\{J[z 2 i-l,Z2i\). 
i=l 



(3-52) 



(3-53) 



These images of the different y-sheets give a partition of the complex plane by definition 
and we can define the boundaries Tj = |J iy separating Q and f2j for any i ^ as 
follows: 

T+ = f]{lmz > 0} = £o+(N-i, z 2i \) = &_(N-i, z 2i \) (3-54) 

and 

r~ = Ti f~~]{Imz < 0} = £o-([z 2i -i, z 2i \) = ^+([z 2 j-i, z 2i \). (3-55) 
We are looking for a solution under the form of a Lax matrix 

/ MM*)) MZiW) ■■■ MM*)) \ 



M(x) 



(3-56) 



where is a Baker- Akhiezer type function analytic on C\ (J^ Tj which must satisfy the 
jump conditions: 



and 



These functions are also constrained by their asymptotic behavior: 

Vi, j = 0, . . . , k, 0i(oj) = 0;(£j(oo)) := ^ 

where one conventionally notes a := oo. 
One also notes 

Tif^R = {p2i-l,P2i} = {€i(z2i-l),Zi(z2i)}. 

One has the solution: 



(3-57) 
(3-58) 

(3-59) 



(3-60) 



MO 



(3-61) 



\ i=i 
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and 



(f>iix) = Ci i ^ (3-62) 



\ i=l 



for i = 1, . . . , k, where the q's are constants which can be fixed by the asymptotic 
behavior of 0j's and the square roots have branch cuts along [J IV 
Indeed, the algebraic equation defining the branch points reads: 

k k k 

- p*-o(£ - a*) = - °*) 2 -E^n^- %) 2 ( 3 - 63 ) 

which gives 



i.e. 



N 

<j>i(a,i) = Ci\/ (3-64) 

Tli 

AT ( ^ 3 " 65 ' ) 

3.2.4 Parametrix at the edge points 

Now the solution of the RH problem can be build away from the branch points thanks 
to the preceding section, let us have a closer look at the branch points Zj. Since the 
behavior near the branch points is the typical square root and one can use the usual 
method (see for example [3]). 

The problem is to find a parametrix in the neighborhood of the branch points, i.e. 
a solution to the RH problem: 

2k 

• P(z) is analytic in ^jD(zj, r)\(R U V) where D(zi, r) is a disc centered at Zi with 

i=l 

small radius r and T is the set lenses defined used for the second transformation; 

• P has jumps: 



r )f| (U D( ^' r) ) ' p+{z) = p -W js W> (3 - 66) 



When n — ■> oo, the boundary conditions for P are: 

Vz = 1, . . . , 2k , Wz G dD{z h r) , P(z) = (l + o(—\ ) M(z) uniformly. 



N 



(3-67) 
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Figure 3: Decomposition of the disc D(z 2i , r) into four sectors. 



Let us now consider a particular branch point ^fl When z — > z 2i , lemma 1331 states 
that: 

Mz) = P2i{z ~ Z2i) + 2j f(z - Z 2i )l +0{{Z- Z 2i f) , v 

A,(*) = p 2 ;(* - z 2l ) - 2ft - 2 2i )f + O ((z - 2 2J ) 2 ) 1 " J 



Ao(^) - A^) = ^(z - z 2l )i + oUz- z 2l )i) (3-69) 



i.e. 

L 

A (z) - Ai(z) = - 

2 

which means that the function (3i(z) = \j(Xo(z) — Ai(z))] 3 is analytic in a neighbor- 
ly 

hood of z 2i with a strictly positive derivative fi[{z 2 i) = p 2i > 0. 

Thus Pi is a conformal map from D(z 2 i,r) to a neighborhood of the origin for r 
small enough. 

Let us now cut the image of D(z 2i ,r) by a a special choice of T representing the 
cuts of the lens based in z 2i : 

A(rf]D(^,r)) = {^|arg(z) = . (3-70) 

It gives a partition of the disc D(z 2i ,r) into four regions labeled I, II, III and IV (see 
figure [3]) following: 

• < arg(/5(z)) < y in region I; 

• ^ < arg(/?(z)) < 7T in region II; 

• —tv < axg(P(z)) < — y in region III; 

• < arg(/3(,z)) < in region IV. 



3 We could do exactly the same for a branch point Z2i-i with the same function /?,. 
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We can now transform the RH for P to a model one by the transformation: 

P{z)Tnr,{z) in regions I and IV. 



P( z ) := i * ppy 

1 P(z) in regions II and III. 



(3-71) 



where one has defined: 



Tp p (z) 



-i k 



'3j 



(3-72) 



J a,0=O 



Thus, P is solution of the same RH problem as P except that its jumps are now given 
by 



VzG(lU 

where jp(z) is given by: 

• for z on [z 2i - r, z 2 i [: 

jp = [(1 — 5 a o — S a i)5 a /3 + 5 a o5f3i — fiai5f3o\a,i3=o 

• for z on the upper and lower sides of the lens in D(z2i, r): 



(3-73) 



(3-74) 



for z on ]zzi, z 2i + H: 



jp= [S a p + S a oS pi e N ^-^] k a 



/3=0 • 



(3-75) 



(3-76) 



One sees that these jumps are equal to the identity except for a 2 x 2 block cor- 
responding to the (0,0), (0, i), (z,0) and (0,0) coefficients. The RH problem reduces 
then to a model problem of size 2x2 which can be solved explicitly as follows (see [3] 
for example): 



(3-77) 




where ensures the boundary conditions 

E N = y/nMMM, 



with 



M = [(1 - S ai (l + i))5 a p - SaoSpi - i5 ai 5po\ k a /3=Q , 



(3-78) 
(3-79) 
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and 



M 



J a,/3=0 



(3-80) 



$(z) is denned with the use of the Airy function Ai(z): 



2tt 



[(1 - <5oo(l - Sfo(z)) - M 1 + ^(^))))*a/3 + S ai S/3oy'o(z) - 5 ai 8 pQ y 2 {z)\ aj3=Q for < arg < , 



,/3=0 



for 



2tt 



< are; < n 



[(1 - 5 a0 (l + y 1 {z)) - 5 ai (l + y' 2 (z))))5 a f3 - 5 ai 5p y'i(z) - S ai 8p y2(z)] a 

[(1 - S a0 (l + y 2 (z)) - S ai (l - y[(z))))S al3 - 5 ai S po y 2 (z) + S ai S m yi(z)} k a/3=Q for - tt < arg < - 
k [(! - *oo(l - yo(z)) - 5 ai (l - y[(z))))5 a fs + 8 ai Spoy'o(z) + S ai 5p yi(z) ]^ =0 for - ^ < arg < 

(3-81) 

where 

'2a.iT; . '2a.iT; . 

Va = e 3 yk(e 3 z) 



2n 



for a = 0, 1, 2 and D is the diagonal matrix 

D{z) = \(1 - 5 a0 (l - e T(*>«-M*))) _ _ e -f (*>(*)-*«(*))) ^ 



a,/3=0 



(3-82) 



(3-83) 



3.2.5 Final transformation 

We now have a parametrix away from the branch points (M(z)) and in their neigh- 
borhood (P(z)). It remains to put them together in a unique function defined on the 
whole complex plane: this is the last transformation. 
Define the matrix 



R{z) = S(z)M(z) 1 for z outside the disks D(zi, r), i = 1, . . . , 2k 
R(z) = S(z)P(z)^ 1 for z inside the disks D(zi, r), i = 1, . . . , 2k. 



One can easily check that: 

R{z) =1 + 



N(\z\ + 1) 



/2k \ 

as N — > oo uniformly for z G C\ | \^jdD(z i: r) J and 



R'(z) = O 



It follows the important relation: 



N(\z\ + l)J ' 



R-\x 2 )R{ Xl ) = I + O (^^) 



(3-84) 



(3-85) 



(3-86) 



(3-87) 



needed to study the behavior of the Kernel in the different regimes. 
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3.3 Mean density of particles 

With the use of the kernel K^, one can state more precisely the role played by the 
spectral curve. It follows from: 

Theorem 3.2 The limiting mean density of eigenvalues of the random matrix 

p(x) := lim (3-88) 

exists, is supported by the real cuts {[221-1, ^2i] }i=i an d is expressed in terms of the 
solution £ of the algebraic equation 

S(x,O = (3-89) 

by 

p(x) = -\ImUx)\- (3-90) 

7T 

It is analytic on the cuts and vanishes like a square root at the edges z%. 
proof: 

Let us consider x e] 221-1 j for any i — 1, . . . , k, assuming that x does not belong 
to any D(z j} r). Thus, Eq. ( 155311 and Eq. fl3^B7j) imply that 

S- 1 (y)S + (x)=I + 0(x-y) (3-91) 

uniformly in iV when y — > x. The expression of the kernel in terms of S(x), Eq. (13-461) . 
gives: 

e N[h(y)-h(x)] 

Kn(x, y) = 2in(x-y) 0l ' s ' i(y) [I + ° (x ~ y)] n isA x ) 



e N[h(y)-h(x)] 



(3 - 92) 



sin (Mm (X 0+ (x) - A + (y))) + 



ix(x - y) 



uniformly in N. Remember that, by definition, 



1 dlmX 0+ (x) 

P( x ) = j • 3 - 93 

7r ax 

Thus, we get 

K N (x, x) = Np(x) + 0(1), (3-94) 

proving the first part of the theorem. For x outside of the cuts (or at a branch point), 
we easily see that 

Kfj(x,x) 

lim NK ' ; = 3-95 

TV^oc N 

because of the ordering of the A functions (see th l3.ll) . □ 
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3.4 Behavior in the Bulk 

Let us first study the asymptotic of the kernel inside the support of eigenvalues. We 
recover the usual sine-kernel: 

k 

Theorem 3.3 For any x G [Jj^i-i) z 2i[ an d any (u, v) G R 2 , one has 

i=i 



1 f> / u r 

hm — ; s K N \ x + ,x + 



iv^oo Np(x) 



Np{x) ' N p{x) ) n(u — v) 



sin TT[U — V) 



where we use the rescaled kernel: 



K N (x,y) = e N ^- h ^K N (x,y). 

proof: 

Consider x &}z 2 i-i, z 2 i[ for any i = 1, . . . , k and define 



X = x + 



u 



Np(x) 



and Y = x + 



Np(x) 



Then, there exists t = x + O (4h such that 

iVIm(Ao + (X)-Ao + (y)) 
From Eq. (j3=52j) . it follows 



7r(u — v)p(t) 

p(x) 



K N (X,Y) 
Np(x) 



(3 - 100) 



sin(Mm(A 0+ (X)-A 0+ (y))) | Q f\ 

ix{u — v) \N 
sin(7r(w - v)) / l\ 

n(u-v) \N J ' 



(3-96) 



(3-97) 



(3-98) 



(3-99) 



□ 



3.5 Behavior at the edge 

Let us now study the behavior of the kernel at the edge of the spectrum, i.e. near any 
branch point of the spectral curve. 
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Theorem 3.4 For any i = 1, ... ,2k and every (u, v) G M. 2 , the rescaled kernel satis- 
fies: 

km ^ Zi + a-, Zi + j 1 



(3-101) 

where pi is a constant characterizing the behavior of p near the branch point z^: 

p{x) = 2i\ x - Zi \k (l + o(l)) (3-102) 

7T 

w/ien x — > Zi. 
proof: 

Let us consider an arbitrary branch point§| z 2 i and define 

w v 
X := z 2 i H 2 an d ^ := z 2i H r (3-103) 

(p 2i A^)3 (p 2i A^)3 

for any (zt, v) g]0,oo[ 2 so that, for N large enough, X and Y lie inside the circle 
D(z2i, r) on the cut. The rescaled kernel can thus be written: 

k n (x, y) = 2 ^_ Y) n 1 MY)s+\Y)S4x)n 2 ^ l (x) (3-104) 

with the notations of Eq. (13-461) and 

S+(X) = R(X)P+(X) = R{X)P + {X) (3-105) 
since we are in the region II. Using the form !3-7Tb f P, one gets: 

S+(X) = R(X)E N (X)$ + (iV§A(X)) D i>+ (X) (3-106) 

where 

(i-l k-i \ 

e JVflm[Ao+(X)] j ^-^ 7e -iViI m [Ao + (X)] ) ^-^ J _ 

It follows 

—L^i^x, y) = -— - — -m- 1 (ivlA(F)) i? Jv 1 (y)^ 1 (y) J R(x)i? JV (x)<i. + 

(3-108) 



6 We consider here a branch point on the right side of a cut but we could do exactly the same 
derivation for the other end point of the cut 221-1- 



28 



where we use the vector of size k + 1: 



A := 



i-l 



k—i 



and 



A := 



•1 ... 1 ... 



i— 1 k—i 

i cT^"o 1 0~^0 



Remember that (3i(z 2 i) = and f3' i {z2i) = /o|j- Thus 

$+ (iVi/^X)) -> $+(«) when N -> oo 



and 



jhn$ + (iV3^(X))^ 



i-l 



k—i 



-i f 



S/o(«) ... ^(m) . . . 



as well as 



lim A*; 1 (iVl/^Y) ) = -1m 



i-l 



k—i 



-y' (u) ... y (u) ... 



with the function y defined by Eq. (13-821) in terms of the Airy function. 
Since one knows the behavior of En, Ay 1 and R when A" — > oo: 



E N (X) 

E N \Y) 
E N \Y)E N {X) 

R-\Y)R(X) 



o im 

(n* 

1 + o (n~* 

I + o(N~i 



(3 - 114) 



we finally prove that 



1 



lim 

(p2iN) 



K N (X,Y) 



yo( u Wo( v ) - y'o( u )yo( v ) 

u — V 



(3-109) 



(3-110) 



(3-111) 



(3-112) 



(3-113) 



(3-115) 



which gives the theorem for u,v > 0. 

For the other signs of u and v, we can derive the result in the same way by using the 
appropriate expression of the kernel as well as the following definition of the h function 
outside of the cut. □ 
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Figure 4: Cuts of the £ functions in the x-plane 



4 Intermediate state 

We consider here one intermediate configuration where two or more segment have 
already merged. 

It corresponds to having at least two non-real complex conjugate branch points. 
4.1 Study of the spectral curve 

We consider an intermediate time in the evolution of the Brownian motions. The 
spectral curve remains the same: 

S(x,y) = (4-1) 

but the value of the coefficients makes its structure different through the position of 
the branch points. 

Let us consider the following case (which is generic). The 2k branch points, solution 
of the equation in y 

i = y — - — 2 ( 4 -2) 

are such that 

• there are 21 real branch points Zi with i — 1, ... ,21 and < I < k; 

• for i — 1, . . . , /, there are < bi pairs of distinct non-real complex conjugated 
branch points (w^ m \n^ m ^j linked by a cut crossing the real axis in a 



unique point rf 1 ^ satisfying: 



\/i , z 2i -i < rf ] < rf ] <...< rf i] < z 2i . (4-3) 
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Remark 4.1 One obviously has the relation: 

i 

k = l + J2bi. (4-4) 
i=i 

and we use a double index notation to sort the filling fraction: 



and the eigenvalues: 



(m) _ (m) . nf m) 1 . 1 I i \ ( 4_7 ) 



Let us now study the properties of the spectral curve as we did in the preceding 
section. We are particularly interested in its sheeted structure. 

A generic point of the complex plane (i.e. not a branch point) has k + 1 pre- 
images on the spectral curve. Let us note these pre- images by £o(x) and $ (x) with 
i = 1, . . . , I and m = 0, . . . , hi. These different pre-images are characterized by their 
asymptotic behaviors when x — > oo: 

C M = a! m) + SjR + O (i) 
From these asymptotic behaviors, one can define the ^-functions by analytical contin- 
uation as follow^ 

i 

• £o{%) can be analytically continued to C\^J[2 2 i-i, z 2 i}; 

i=l 

• for i = 1,...,/ and m = 0, £j"^(a;) can be analytically continued to 
T-" 1 '' 1 U [r l - m ' ) , r™' +1 ] U r^ OT+1 M where we use the notation: 

(0) _(o) (o) 
w \> = w y = r \> := z 2i -i, 

_ ^fti+D _ jbi+i) ._ 

W i — W i — r i ■— Z 2i- 



It implies the sheet structure described in fig{ 

• the 0-sheet is connected to the (i,m)-sheet along the real cut [r\ , r t - m+1 ^] for 
i = 1, . . . ,1 and m = 0, . . . , 6j with the jumps: 

M*) = C ( ?(x) for x G]rJ m) , rf +1) [; (4-9) 



7 For the sake of brevity, I did not include any proof of this sheet structure in the present manuscript. 
Nevertheless, the simplicity of the spectral curve implies a simple, but long and uninteresting, proof 
of these facts. One can, for example, follow the time evolution of the sheet structure from t = 1 to 
t = through the rational parametrization of the algebraic curve. As such a proof does not involve 
any new technical element, it is left to the reader as an exercise. 

8 We label the sheets (except the 0-sheet) by two indices (i,m) corresponding to the associated £ 
function £. . 
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Figure 5: Sheeted structure of the spectral curve for an intermediate step. In this 
example, there are 4 real branch points and 4 pairs of conjugated ones. 



• for % — 1, . . . , I and m — 1, . . . , b{ — 1, the (i, m)-sheet is connected to the (i, m—1)- 
and (i,m + l)-sheets along the imaginary cuts r\ m ^ and r^ m+1 - ) respectively, with 
the jumps: 

&\*)=&r 1 \x) for xer<T\ UW ) 
&\x)=fg ¥1 \x) for xerf +1 ). (4 " 10) 

We finally define the associated A functions as in the preceding section: 
Definition 4.1 Let the function be the primitives of^ m ^ defined by the conditions: 

Ao(*a) = °> ( 4 " n ) 

AS5(*2i-i) = Ao-(z2i-i), (4-12) 

\ < ? i) (z 2i )=\ Q+ (z 2i ) (4-13) 

and 

A M (w ( m+ i )) = A (-+D (w; (™+i) ) (4 _ 14 ) 

/or m = 1, . . . ,bi - 1. 

One can easily prove that they satisfy the following properties: 
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i-1 

Lemma 4.1 Ao «s analytic in C\^J[z2i-i, &n] and A-" 1 " 1 is analytic in C\^J [22.7-1 5 ^2j]U 

2i-i) U I^ m ' ) /or z = 1, . . . , k and m — 1, . . . , Oj. 
TTie jumps and asymptotics of these X-functions are given by: 



When z —>■ 00: 

z 2 , , „ / 1 



A = --mz + / + O(-) (4-15) 



2 U 

and, /or z 7^ : 



(m) , 1 

\< m) = a< m W^ln* + i< m >+o(^) (1-16) 



where Vf 1 's are some constants chosen to satisfy the constraints ( 4-11 ) — ( 4-14 )- 
They follow the jumps: 

1 b J n (m) 

A 0+ (x) - A _(x) = -2m ^ ^2 ~pr> x e l^i-2, z 2i -i[ 

j=i+l m=l 

Ao+(x) — Ao-(x) = —2in, x g] — 00, zi[. 

(4 - 17) 

For z G [1, /] 

(m) 

A^ (z) - AP(x) = 2^, x g] - 00, rl m) [ 

A +(x) - X^y{x) = , A -(x) - Xi+\x) = 0, x e]«a-i, z 2 i[ 

(4 - 18) 

Fori G [1,1- 1] 

A 0+ (x)-AP(x) = 0, xG[r^,rf +1) ] 

(4 - 19) 

and 

X%\x) = X% +1 \x), self** (4-20) 
where I\- m ^ + denotes the part of T • m ' ) lying inside the upper half plane. 

(m) x (ft) 



4.1.1 Study of the sign of Re (\ (m) - Aj 

It is very important to know the sign of Re (^xf™^ — X^j for (i, m) 7^ (J, h) in order to 
study the asymptotic behavior of the jump matrices in the Riemann-Hilbert problem we 
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^(m) 



want to solve. For this purpose, it is convenient to draw the curves where Re ^A- 

Re Let us build them step by step. 

First of all, from the behavior of the A-functions near the branch points, one can 
state the following results: 

• From any branch point z^i-i go three curves at equal angles ^ where Re (Aq) = 



Re one °f them being 



(i) 



From any branch point z 2 i go three curves at equal angles ^ where Re (Aq) 



Re (xf^ j one of them being 



• From any branch point (resp. w\ ) go three curves at equal angles ^ where 
Re ^Aj- m = Re^A^J, one of them being wf^.ioo (resp. w\ m \—ioo ). 

Let us denote r- r ^ + ' ) and r i -™ + ' ) (resp. T-^ ^ and v'fl ■*) the two other curve going 
from wf 71 (resp. w\ m ^), the first one lying to the right hand side of the curve T^' 
(see fig. E]). The curves and ^ (resp. r.™^ and ) intersect in the 

real line in xf^ (resp. xf^). On the analytical continuation (resp. T?^P) 

of r^,^ (resp. T\™ ±S> ) beyond xf^ (resp. x\^), one has Re ^Ao — A-" 1-1 ^ = 

(resp. Re ^Ao — A- m ' ) j = 0) because of the cut on the real line. 

We have investigated the curves Re ^A- m ' ) — A^J near the branch points, let us 
now look at their behavior away from them. Many different configuration can occur, 
depending on the relative positions of the xf^'s, the x^'s, the r^'s and the real 
branch points Z{. In order to limit the size of this paper, we restrict our study to 
the cases where rj < xf"^ < rf 1 +1 ' and r\ m ' < < rl m+1 \ Nevertheless, the 

RH problem corresponding to the other possible configurations can be solved with the 
same method and the same transformations even if we do not give a proof here (the 
study if the signs of Re (xf 11 ^ — A^ j is a little bit more subtle but does not present 
any particular difficulty). 

Note also that the curves Re ( X^A = Re ( X<; m) ) and Re ( X^ n+1) ) = Re ( X {m) 



cross in two points symmetric with respect to the real axis. On the line d\ joining 
these two points, one has Re [X^™ 1 ^ = Re fAj m+1 M. 

With this assumption, between and r[ m+1 \ there are two different possible 
configurations: 

• In the first case: 

-i m) <*S<*S +1) <^ +1) - (4-21) 
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Figure 6: Zoom near the vertical cut T^. The thick solid lines are the curves 
where Re yX\ m ^ — A^ m_1 ^J = 0. On the dashed lines (resp. F^), one nas 

Re (a - A^" 1 *) = (resp. Re (a^ - A ) = 0). 



The respective weights of Re (A? 1-1 (a;)), Re(Af l (x)), Re(A™ +1 (x)) and 
Re ^AQ°^(:r)J for x between and r\ m ~^ are depicted in figJ3 

• In the second case: 

r M < x (rn + l) < x 0») < r ( m+ l)_ (4 _ 22) 

The respective weights of Re (Af -1 ^)), Re(A™(x)), Re (A™ +1 (x)) and 
Re ^Ao^(x) j for x between and r^ m_1 ^ are depicted in figJBJ 

In both configuration one has the following result: 

Lemma 4.2 For x between and rj 7r 1 : 

\f \x) > A (x) (4-23) 

for (j, /i) ^ {i,m— 1), (i, m), (z, m + 1). 

4.2 Riemann-Hilbert analysis 

We now solve the Riemann Hilbert problem presented in the preceding section for 
this particular value of the external matrix eigenvalues by performing a sequence of 
transformation leading to a simple problem when iV — ► oo. 
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Figure 7: Study of the sign of Re (xf ) - for the first case. For readability, 

m > m+l > m-l > stands for Re (A™(x)) > Re (A™ +1 (x)) > Re (A™ _1 (:r)) > 
Re(A°(x)). 



For convenience, we use a double index notation for the matrices in this section: 
considering an arbitrary k + 1 x k + 1 matrix A, we note its matrix elements as follows: 



[A\{i,m);(i' ,m') '■ — [A\^i , ; j. , 



J E;=i 6 J +m+i,E;=i 1 b J+ m'+i' ( 4_24 ) 

for i = 1, . . . , / and m = 0, . . . , hi. The first line (or column) is denoted by the pair 
(0,0). This notation coincides with the notation used to sort the branch points and 
filling fraction as well as the t}™^ and A^ functions: the lines are sorted in i groups 
of bi elements corresponding to the i real cuts \z2i-\, Z2i] shared by hi + 1 sheets. 

With these notations, let us remember the RH problem to solve. We are looking 
foraA; + lxA; + l matrix Y(x) satisfying the constraints: 

• Y is analytic on C \ 1; 

• for x G K, one has the jumps 



Y+(x) = Y-(x) 



( 1 




V o 



^i 0) (*) 






1 











1 











1 
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(4-25) 



m-l=m m-l-m+1 

m>m-l>m+l>0 



m=m+l 



m+l=0- 



m-l>m>0>m+ 



m>m+l>m-l>0 




m-l=C 



Figure 8: Study of the sign of Re — A^ 1 ^ for the second case. For readability, 

m > m + l > m-l > stands for Re (A™(x)) > Re (A™ +1 (:r)) > Re (A™ _1 (x)) > 
Re(A°(x)). 



where Y + (x) and Y_(x) denote respectively the limit of Y(z) when z 
the upper and lower half planes and ujf^^x) := e N ( 2 ai x ). 



when x — > oo, one has the asymptotic behavior: 

/ z N ( 



Y(x) 



i + o ( - 

x 



Z~ n i 



(0) 



z~ n i 

V o o o 



(1) 



o \ 






Ah) 



J 



4.2.1 First transformation 



x from 



(4-26) 



The first transformation is the same as in the large time case when the branch points 
are all real. 
We define: 

T(x) := n ltYT Y(x)n 2tYT (4-27) 

l bi 

for x e C\R(J |J K (m) ,tZjS m) ], with 



j=l m=l 



^^diag^^-^,...^-^) 



(4-28) 
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and 



il 2Y T ■= diag eV 2 J ,e \ 



,...,e 



(4-29) 



Then one can check that it is solution of the Riemann-Hilbert problem: 

l h 

• T is analytic on x G C\M^J 



(m) — (m)i 



i=l m=l 



T has the jumps 



T + (x) = T_(x) Jr (x) 



where 



/ e iV(Ao+(z)-Ao-0z)) e N(\[°l(x)-X .(x)) e iV(AW(x)-A _(x)) 





e iV(A«(*)-A«(x)) 















(4-30) 



e JV(A#>(x}-Ao-(x)) \ 






^(a^'^-a'^W) 



(4-31) 



for x G R and 



5, 



for x G [u; l - m+1 \ u;j m+1 -'] with i — 1, . . . ,1 and m = 0, . . . 6j — 



/ ArfA^-A^) 
1 + {d{j,h);(i,m) + <%/i);(i,m+l)) ( e V J+ — 1 

(4-32) 



when x — ► oo, one has the asymptotic behavior 

1 



T(x) =1 + 



x 



(4-33) 



Note that one can factorize the jump matrix Jt{x) depending on the cut on which 
x lies: 

For x G [rj OT \ r^ m ] with i = 1, . . . , I and m = 0, . . . , bi, all the diagonal terms 
of Jt reduce to 1 except the first one and the term [(i,m); (i,m)] whereas the term 



9 It simply means that the jump matrix is the identity for x <G [u>i > w i ] except for the 

(«,m)'th and (i,m + l)'th diagonal terms which are equal to e v i+ *~ J and e V i+ i_ / 
respectively. 
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[(0,0); (i, m)\ (i.e. the (z,m)'th term of the first line) also goes to 10: 

/ e 7V(A 0+ -Ao_) p iV(A^-Ao-) 










e K !+ 
1 








(m-l)_ A \ 


1 


Ar( - A (m + 1)_ A , 


N(\<; bl) 














1 














e N{x^-xt ] ) 














1 














1 



(4-34) 



For x e M\ Uit^i-i, ^2i]j all the diagonal terms of jx reduce to 1 and it takes the 
form 



Jt 



( l e N(\[°_l(x)~\ -(x)) e 7V(A^(x)-A -(x)) 

1 

1 



\ o o 

The kernel can be written: 







^(a&'Vj-Ao-W) \ 




1 / 



(4-35) 



K N (x,y) 



4 4 



fi liT ( 2 /)T- 1 ( 2/ )T + (x)^ 2iT (x) 



where 



and 



2in(x — y) 
n liT (x) = [0 ... e^^] 

fi 2T ( x ) = [ e -^o+W o ... 0]. 



(4-36) 

(4-37) 
(4-38) 



4.2.2 Second transformation, global opening of lenses 



We now proceed to the second transformation allowing to cancel the exponentially 
increasing coefficients of the jump matrix jj<. This transformation was introduced in 
pE], for the particular case of two imaginary branch points i = 1 and b± — 1. We 
generalize this procedure using intensively the study of the spectral curve presented in 
the preceding section. 

i 

For this purpose, we define a contour £ = [^JSj composed of I connected Compo- 
te 

nents build as follows (see figures [3, [8] and M) : 
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We show here the form of jr for m ^ 0, hi. These two cases take exactly the same kind of form. 
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Figure 9: Example of contour E used for the global opening of lenses. Here, h = 2 and 
we are in the first case where xf^ R < xf\. 

• Ej is a closed curved encircling the cut [z 2 i-i, z 2 ^\ cutting the real axis in X2i-i 
and x 2 i satisfying: 

X2i-2 < x 2i -i < z 2i -i < z 2i < x 2i ; (4-39) 

• Ej intersects the line rf^ + iR in wf^ and wf^ for m = 1, . . . , bf, 

• in a neighborhood of wf^ and w\ m \ Ej is the analytic continuation of the curves 
Re (A l (m) )=Re(Af- 1) ); 

• in the half-plane to the left (resp. to the right) of rf^ + iM., E lies in the region 
Re (Af > Re (aJ™ } ) (resp. Re (a^ < Re (\ (m) ) ) for any m = 1, . . . , ft*. 

Let us now define the second transformation T — > [/. 

Definition 4.2 Let f/(x) fre £/ie A; + 1 x A; + 1 matrix defined by U(x) = T(x)Jut(x) 
where Jut is given by 

• For x outside E ; Jjjt(x) = I; 

• For x inside Ej and Re(x) < J"j(l): 

r 7 / M r r r v(A (1) (a;)-A (()) (a;) N ) ,„x 

i J ^(»)J(,-,h);( J -' )V) ■=0(j,h);U',h') -5{j,h);{i,Q)5{j',h');{i,l)e ^< >\ (4-40) 
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For x inside £, between T ■ m ' ) and H m+1 ^ : 



[ J UT{x)} {jA y Xfih/) 



n( A (m_1) (x)-A (m) (x) 
°(j,h);(j',h') — °(j,h);(i,m)0(j' ,h');(i,m-l) e ^ ' 

— 0(j,h);(j,m)tf(j',h');(i > m+l) e i 



(4-41) 



For x inside Sj and Re(x) > ri(bi): 



[JuT{x)] {jmjl)hl) := <*(,-,&);(,•',*»') - 5( J> );(i,6 i )%',h');(i,6i-l) e 



(4-42) 



The matrix [/ (x) is then solution of the Riemann-Hilbert problem: 



U{x) is analytic in C\ ( R \J S|Jr 



(m) 



1,7X1 



when x — ► oo, one has the asymptotic behavior: 



U{x) 



! + ()[- 

x 



(4-43) 



One has the jumps U + (x) = U-(x)ju(x) with the jump matrix given by: 

+ for x e]r , ?4 m+1 ^ [, with z = 1, . . . , / and m — 0, . . . , 6j, the jump of the matrix 
{7 is the same as jump matrix of Eq. (14-341) except the (i,m — l)'th and the 
(i, m + l)'th terms of the first line which are taken equal to 0: 



1 



N(*i+\*)-*i'"'(.°>)) 




e JV(A 1 ( 'j i) (x)-A _( :c )) ^ 









(4-44) 



+ for x G [z 2 *-i,7i ]: 



/ e 7V(A 0+ (x)-A -(x)) 






1 e"™+ 

e 7V(A(°)(x)-A(°)) q o 

1 



JV(A$(aO-Ao-(<n)) g^C^ (*)-*>- (*)) \ 






1 

(4-45) 
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+ for x G [rf i \z 2 i} ) ju is given by: 



( e N(A 0+ (*)-A _(*)) e N(AW(*)-A _(*)) 
1 



+ for x G R\ Ui [x 2 i-i,X2i] 








e «(A< t [) 1+ (x)-A _(x)) 






1 



N(A <! ;' ) (x)-A _(x)) \ 



(4-46) 



/ 1 e N(X^(x)-\ {) -(x)) e iV(A^(a!)-Ao-(x)) 






V o" 



1 







1 



e JV(A^ ) (x)-A _( ;c )) \_ 






/ 



(4-47) 

+ for x G [x 2 i-i, z 2 i-i], ju is the same matrix as for x G M\^J[x 2i _i, a^i] with the 

i 

term (i, 1) of the first line taken equal to 0. 

+ for x G [z2i, X2i\ ju is the same matrix as for x G R\^J[x2i-i, X2i] with the term 

i 

(i, bi — 1) of the first line taken equal to 0. 
+ for x G r (m) 



\ju(v)\{j,h);{j',h') 



1 - <%fc);(i,m) + <%fc);(i,m-l) e 



A? 



( A ( T -D_ A ( r D) _ x 



+^0»;(i,m)%'^');(i,m-l) ~ %^);(i,m-l)%'^');(i,m) 



— S{j,h);(i,m-l)^{j',h');(i,m-2)e 



N 



^ A (™-2)_ A (m-l)~J 



+<J(j,fc);(i,m) 



<%',/i');(i,m+l) e 



N 



V ; - 0(j',h');(i,m-2)e 



m — 2) ^(^) 



(4 - 48) 

It means that the jump matrix reduces to the identity outside the 2x4 block: 



(i, m — 1) 
(i,m) 



(i, m — 2) 



3 Af(A( m - 2) -A<7') 



(i,m— 1) (i, m) (i, ra + 1) 



-1 

r 







Ar(\ (m+1) -A< m) ) 



(4-49) 



+ for x G Ej between x 2 i-i + iM. and + iR 

\ju( x )h,W,h') : = 5 V,hW,h') + 5 U,h);(i,oAj',h');(i,i)e N ^ {x) ~ K (a °) ; (4-50) 
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+ for x G Ej between and r j - m+1 ' ) : 

r • i w r r ^ Jvf^ m - 1) (s)-A t ( m) ( a! )') 

:= °(j,h);(f,h') + d (j,hy,(i,m)0(j',h');(i,m-l)e V ■ 7 

. <■ , w(A^ m+1 ^(z)— A^ m ' (x)) 

+C , 0',/i);(j,m)"(j',/i');(*." 1 +l) e ' 

(4-51) 

+ for x G Ej between + iR and X2« + iR 

N(x {bi ' 1) (x)-X {bi) (x)) 

\ju( x )\(j,h);(j>,h') : = S U,h);(3',h') + S (j,h)l«b i ) S (j',h');{i,b i -l) e V * * A (4-52) 

Using the study of the sign of Re(A- m ' ) ) — Re(A^ ) performed in the beginning of this 
section, let us check that these jump matrices have no exponentially increasing entries 

when TV — ► oo. For x G E and x G 



°\^J[02i-i) Z2i]\ j the jump matrices converge 



to the identity matrix as N — > oo. For x G v\ , the jump matrix 14-481 reduces to the 
identity except the block 0321 which converges to the block 





(i, m — 2) 


(i, m — 1) 


(i, m) 


(i,m + 1) 


(z, m — 1) 








-1 





(i, m) 





1 









(4-53) 



exactly as it was encountered in the simpler case studied by p|]. 

Finally, all the non-constant off-diagonal entries of the remaining jump matrices 
!4-44[ l4^4"5l and |4"-46I converge to as iV — > oo. The non- constant diagonal entries have 
modulus 1 and are rapidly oscillating as N — ► oo. The third transformation aims to 
turn them to exponentially decreasing terms. 

4.2.3 Third transformation 

The third transformation aims to cancel the oscillating parts of the jump matrices by 
opening lenses on [z 2 i-i, z&] as in the large time case. 

Remark hat one can further factorize the matrix ju{x) for x G [r\ , r- m+1 ']: 

ju (x) = jus (x)js (x)jus (x) (4-54) 
where, for x G [r^, 7$ ], one has defined the k + 1 x k + 1 matrices: 

l(js)](j >h y,(j>,h') := 0- ~ &(j,h);{i,m) ~ %fe);(0,0))%ft);(i',ft') + 

+8(j,h);(0,0)fi(j' ,h');(i,m) ~ %,h);(i,m)<%',h');(0,0) + 
V- , AT^-Ao- 



+S(j, h y,(o,o) }^ s (j',h');(j",h") e 

U",h") jt (0,0), (i,m), 
(i, m — 1), (i, m + 1) 
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Figure 10: Example of lens for the third transformation. 



-6, 



(j,h);(i,m) 



^2 S U',h');U",h") e 



Y(A^-A 0+ 



(4 - 55) 



U",h") ? (0,0), (i,m), 
(z, m — 1), (i, m + 1) 



and 



(j,h);U',h>) 



:=6, 



(j,h);(j',h') + 8(j,h);(i,m)0(j',h');{0,0) e 



,AT(Ao-A< m) )- 



X I r r W(A -A (m) ) + 



(4-56) 



(4-57) 



(j,h);(j',h') 

As in the large time case, we define a set of I lenses Lj with vertices z<n and ^i-i 
(see fig JTOl) which do not intersect the contour S. Using these lenses, we define the 
change of variables: 

S(x) := U(x) (4-58) 
outside of the lenses and, for x inside the lenses: 



S{x) 



U(x)jjjg\x) in the upper lens region 

T (do 

'us 



U (x) jjjg Wn \x) in the lower lens region 



with 



J, 



(up) 



'■= S (j,h);(j',h') - %ft);(*,mAj',/ l ');(o,o)e Ar(Ao x * )} 



J, 



(down) 
US 



■= hhhW,h') + $(j,h);(i,m)5(j>,h>);(0,0)e N(X ° ^ >] 



(3,h);(3',h>) 

for x between r„- m ' ) and r„- m+1 ' ) 



J, 



(up) 
US 



(j,h);(j',h') 



:=6, 



(j,h);(j',h') ~ S(j,h):(ifi)°(j',h'y,(0fi) e 



N(\ Q -\f ] ) 



J, 



(down) 
US 



U,h);U',h>) 



:= 5, 



(j,h);(j',h') + $(j,h);(i,0)°(j',h>);(0,0) e 



5ai h'\ 



^(Ao-Af 5 ) 



(4-59) 



(4-60) 



(4-61) 
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for x between z 2 i-\ and and, 



J, 



(up) 



${j,h)-,(j',h>) - $u,hy,(iA) 6 ti',h'y,m eN(Xo x<i ' )} 



J, 



(down) 
US 



. w ., .„ := S m;(f,h') + S(j,h);(.iMAj',h');m e 

U,h);(y,h') 



N(X -4 bi) ) 



(4-62) 



for x between rf 1 ^ and z 2i . 

We finally get the jump matrices for S(x) out of the discontinuities of U (x) thanks 
to the factorization of the jump matrices for U on the real axis. One has the jumps 
S+(x) = S-(x)js(x) with the jump matrix given by: 

„(m) (m+1) 



for x e]r,- , r\ m [, with % = 1, . . . , I and m = 0, . . . , bf. 



[Us)} 



U,h);0',h') 



(1 - S(j,h);(i,m) ~ S(j,h);(0,0))8(j,h);(j',h') + <%fc);(0,0)%',fc'); 



);(i,m) ~ S(j,h);{i,m)S(j', 



+ S (j,h);(0,0) 



>(j,h);(i,m) 



E 

t (o,o; 

E 



(3",h") # (0,0), (i,m), 
(i, m — 1), (i,m + 1) 



%',/>');(i"^") e 
%',h');0",ft") e 



AT A 



(4 - 63) 



U",h") (0,0), (i,m), 
(1,171 — 1), (i, m + 1) 



for x G R\Uj[x 2 i-i,a;2j 



/ 1 e JV(A^(x)-A () -(x)) e JV(AW(a!)-Ao-(x)) 

1 

1 



V o 











e N(\^ ) (x)~X -(x)) \ 






/ 



(4-64) 



• for x G [x 2 i-i, z 2 i-i], js is the same matrix as for x G K\ [x 2 i-i, x 2i \ with the 
term (i, 1) of the first line taken equal to 0. 

• for x G [z 2i , x 2 i] js is the same matrix as for x G R\ Uj [x2i-i, x 2 i] with the term 
(i, bi — 1) of the first line taken equal to 0. 

• For x G Tf 1 ^ between x\"^ p and and between wf^ and xf^ own 



\Js(x)] 



U,h);(j',h>) 



\hh);(j'X 



1 - $(j,h);(i,m) + $(j,h);(i,m-l) e 



N 



-1) .(m-l)^ 



+%^);(j,m)%',/i');(i,^-l) ~ %,ft);(j,m-l)%',h');(i,m) 
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~ 5{j,h);{i,m-X)5{j' ,h');{i,m-2)^ 



JV( A (m_2) -A (m_1)N 



(j,h);(i,m) 



5(j',h');(i,m+l)£ 



(m+l) ->(m)\ 



ATI A; —A 



3(j',h'y,(i,m-2)e 



JV A 



(m~2) 



(4 - 65) 

It means that the jump matrix reduces to the identity outside the 2x4 block: 
(i,m — 2) (i,m— 1) (z,m) (z,m + l) 



(?j m — 1) 

(i, m) 

-,(m) 



Af(A (m - 2, -A (m - 1) ) AT(A^ _1) -A^ m_1) 
7vfA (m - 2) -A ( ™') 



-1 






jvfA (m+1) -A (m)N 



(4-66) 



For x G Tj 1 " 1 ' 1 between r} and 



™( m ). 



(j,h);(j',h f ) 



6 



{j,h);(j',h>) 



1-6, 



(j,h);(i,m) + fyj»;(t,m-l) I e 



r r JvfAo-A^" 1 )) 

" d 0'.'»);(».»w-l) <) (j / .'i / )i(0.0) e 



+O(j,/i);(l,m)0(;7' ) /i 



(3,/i);(i,m)O0' / »ft / );(*.»n-l) ~ °a,h);(i,m-l)%',>y);(*>w) 



%,h);(i,m-l)0(j',h');(i,m-2) e 



Ni A^ m " 2) -A (m_1)N 



+5 



(j,h);(i,m) 



(4 - 67) 



./vYa?" 1 " 1 " 1 - 1 — A? m ^ TA/Y 

%',/i');(i,m+i)e v 1 y — <5(j/ ) / l ') ; ( i>m _2)e v 



A (m-2) 



i.e. this is the same jump matrix as 14-651 except the term (0, 0) of the (i, m — l)'st 

jvY An— A^ m_1 'l 

line which is equal to-e ' ^ / instead of 0. 



For x G r. m ' ) between £-™]„, n and r- m ^: 



(j,h);(j',h') 



5(4 



1-6. 



N{X -\) 



(j,h);(i,m) "I - 0(j»;(i,m-l) I ' 

(m-l)\ 



+%,/i);(i,m-l)%',/i');(0,0)e 
Jr 5(j,h)-(i,m)5(f ,h');(i,m-l) ~ 5(j,h);(i,m-l)5(j' ,h');(i,m) 

— 5(j t h);(i,m-l)5(j>,h');(i,m-2)e 



Nl A (m_2) -A (m_1) 



+5, 



(j,h);(i,m) 



(4 - 68) 



, jv(A< m+1) -A (m) ) r n() 

0(j>,h');(i,m+l)e V 7 — 0(j>,h');{i,m-2) e V 



(m-2) 



i.e. this is the same jump matrix as 14-651 except the term (0, 0) of the (z, m — l)'st 

jV('Ao-A (m_1 ^ 

line which is equal toe ' i_ / instead of 0. 
for x G Sj between a^j-i + £R and rj + z'M 

Vs( x )\(j,h);(j>,h>) : = °U,h);U',h') + %fc);(i,Q)<%',ft');(i,i) e v 1 1 y ; (4-69) 
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• for x G Ej between and r| m+1 ^: 

:= d U,h);U',h') + d O>);(i,m) d Cj',ft');(i,m-i) e v ' ; 

. . , iv(A< m+1) (*)-A< m) (s)) 

+ d (j.'»)i(*.n») d Cj'.h')i(i.n»+l) e V ^ 

(4 - 70) 

• for x G Sj between rf^ + and a^i + iM. 

• for x G U between r\ (m) and I\ (m+1) 

N(X A^ m ^ ") 

[ ,7s ](3'.'»);(3''.h') := % /l );0''.' 1 ') + %M;(^)%'^');(o,o)e °~ 1 (4-72) 



for x E Li between zn-\ and T- 1 ^ 



[ J s]{j,h);(j',h') '■= S (j,h);(j',h>) + %,/ l );(i,o)^(i',fe');(o,o)e JV(Ao A * )} (4-73) 

• for x G Lj between Tf^ and 2 2 i 

[^lo'.^Ki'./i') := 5 0,h);U',h') + 5 (i,/ l );(i,b ! )%',/i');(o,o)e JV(Ao_Al ' } (4-74) 

Now, when N — > oo all the non-constant terms of the jump matrices tend to 0. We 
are then left with a model Riemann-Hilbert problem reduced locally to a 2 x 2 matrix 
problem similar to the one studied in the preceding section. We can thus use the same 
parametrix. 

4.2.4 Model Riemann-Hilbert problem 

We proceed exactly as in the preceding section (large time case) and consider the 
following RH problem (this is the limiting problem obtained when N — > oo). 

Definition 4.3 M is solution to the following RH problem: 

• M : C\ (Uifza-i, *2i] U,, m rS m) ) -> C fe+lxfe+1 arm/^c 

• M(x) = h+ixk+i + O (±) when x -> oo 

• M(x) satisfies the jumps M + (x) = M_(x)j M (x) with: 

(4-75) 
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for x E]rf n \rf n+1 ' [ and 

frM( X )]( ? >),G>',ft / ) = S (J,h),(J',h') (l - <%/0,(»,m-l) - %/i),(i,m)) + %/ l ),(i,m-l)£o',fc'),(»,m) 

-%',h'),(i,m-l)%h),(i,m) 

(4 - 76) 



/or x G 1 1 ■ 

We proceed as in the large time case, following [4J, to solve this problem: we build 
a solution under the form of a Lax matrix associated to the spectral curve. For brevity 
we will not enter the details of the resolution but only give the result, the proof is 
similar to the building of the parametrix away from the branch points in the large time 
case. 

The matrix M takes the following form: 

(h) (Ah 1 ), 



mx)) m ,vM = M) n >\?f>{x)) (4-77) 
where the functions M- are Baker-Akhiezer functions defined by: 



Mf (x) = ) ' (4-78) 

y/R{x) 



and 



where 



H(i.m)^(i.h) ( x a 



where 



and 



(x) = c ' — — (4 " 79) 

<!<"" := ?f'(oo), (4-80) 

I 

R(x) := fj {x - P 2i-i) (x - Vii) J] - ir ] ) - i n) ) (4-81) 

i=l i,m 

<lt ] ■= Co (^! m) ) and ■= £o (^ (m) ) (4-82) 



(m) 

c! m) :=-n/V- (4 " 83) 

We now look at the behavior of the solution in the neighborhood of the branch 
points by approximating it by a local paramerix build from the Airy function. 

Near the real branch points Zj, we use the local parametrices P built in the preceding 
section (see Eq. (13-771) ). Around the imaginary branch point, we have to build a slightly 
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different parametrix. For this purpose, we can directly use the parametrix built in 
section 7 of []. Let us remind here this result for the parametrix around w^ m \ 

Around the branch point let us define a local parametrix P{x) for x inside 

a disc T>^ m \r) of small radius r centered in w^ m \ Studying the behavior of \\ m \x) — 
\\ m ~ 1 \x) inside this disc, one can define a conformal map f(x) by: 



/(*) == 



( \ ( m ) i \ \ ( rn ' 
-(A) \x)-\\ 



(4-84) 



One can remark that the neighborhood T>\ m \r) is mapped to the complex plane in 
such a way that: 



arg(/(x)) = < 

3 



2tt 
3 



for x e r\ (m) 

for x G Ej in the left half-plane 
for x G Sj in the right half-plane 



(4-85) 



With these tools, one defines the parametrix P(x) in T>\ m \r) by: 

P(x) := E(x)$ (N*f(xj) D{x) 
with D{x) is a diagonal matrix defined by 



(4-86) 



i D ( X )](j,h);(j',h') : = 5 (.j,h);(j',h>) 



1 + <%fc),(i,m-l) e 



N ( Am 



- lj + <5(j»,(i,m) 

(4-87) 

for x G Pj- m ' ) (r), $(x) is the identity matrix except the 2x2 block corresponding to the 
intersection of the (i,m — l)'th and (i,m)'th lines and column which takes the value 



for < arg (f(x)) < 



2tt 





(i, m — 


1) 


(i,m) 


(i, m — 1) 


-y' 2 




y'o 


(i,m) 


-2/2 




?/o 




(i, m — 


1) 


(i, m) 


(i, m — 1) 


-y[ 




y'o 


(i, m) 


-2/i 




2/o 



(4-88) 



(4-89) 



for > arg (f(x)) > and 





(i, m — 1) 


(i, m) 


(i, m — 1) 


"2/2 


y[ 


(i,m) 


-2/2 


?/i 



(4-90) 
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for ^ < arg (/(#)) < ^ with the functions y a defined in Eq. (13-821) . The matrix E 
is a prefactor used to ensure that the local parametrix can be linked properly to the 
global parametrix on the boundary of T>^ (r) : 

E(x) := M(x)L(x)- 1 where L(x) : = —=L(x)L(x) (4-91) 

with L the diagonal matrix 
H x ) := 5 m . {j , thl) \l + <J a ,ft);(i, m _i) (Nvf(x)i - l) +6 {jmLm) (N-vf(x)-i 

1 J U,h);(y,h') L \ / \ 

(4-92) 

and L the identity matrix except for the 2x2 block corresponding to the intersection 
of the (i,m — l)'th and (i, m)'th lines and columns: 





(z, m — 1) 


(z, m) 


(z, m — 1) 


z 


-1 


(z, m) 


z 


1 



(4-93) 



4.2.5 Final transformation 



We now have all the ingredients in hand to perform the final transformation and build 
a solution to our Riemann-Hilbert probem. Let us define: 

R(x) := S'(a;)M~ 1 (x) outside of the discs surrounding the branch points , . Q 
R(x) := S(x)P(x)~ 1 inside of the discs surrounding the branch points 

where P(x) is the local parametrix Eq. (13-771) if one considers a real branch point or 
Eq. (14-861) if one considers a complex branch point. 

As iV — > oo, it is solution to the following RH problem: 

• R(x) is analytic inside the discs surrounding the branch and in 

c\ iRU£U^U r * M ) outside; 

• its jumps outside the discs are of the form: 

R+ = R - (I + O (e~ aN )) (4-95) 

with a > and its jumps on the circles surrounding the branch points are such 
that: 

R+ = R-(l + o(±y). (4-96) 

• as x — » oo: 

R(x) =1 + f^Y (4-97) 
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As in the large time case this implies the behavior of R: 



R(x) =1 + 



( 



N(\x\ + 1) 



1 



) 



(4-98) 



as N — > oo (see jl]). 

4.3 Asymptotics of the kernel and universality 

Using this final transformation as well as the properties of the parametrices, one easily 
shows theorems 11.11 11.21 and 11.31 by mimicking the proof of the large time case. 

5 Perspectives 

We used the Riemann-Hilbert approach to study the asymptotic behavior of N non- 
intersecting Brownian bridges and show the universality of their correlation functions 
as N — > oo. Nevertheless, we did not study the critical configurations corresponding to 
a spectral curve where two real branch points merge. One expects to find a universal 
behavior given by the Pearcy kernel as it was found in the particular case studied by [5] . 
This critical case may be addressed by the same techniques and is left to a forthcoming 
work. However, as the length of these notes witnesses, the RH approach is long and 
requires a global information on the spectral curve to finally give local results. It would 
be interesting to obtain these results in a more direct way using, for example, a direct 
steepest-analysis as performed by Adler and Van Moerbeke [lj. 
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